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Adequate description of electro and gamma nuclear physics is of utmost importance in studies of electron 
beam-dumps and intense electron beam accelerators. I also is mandatory to describe neutron backgrounds and 
activation in linear colliders. This physics was elaborated in Geant4 over the last year, and now entered into 
the stage of practical application. 

In the Geant4 Photo-nuclear data base there are at present about 50 nuclei for which the Photo-nuclear 
absorption cross sections have been measured. Of these, data on 14 nuclei are used to parametrize the gamma 
nuclear reaction cross-section The resulting cross section is a complex, factorized function of A and e = log{E 1 ), 
where Ej is the energy of the incident photon. 

Electro-nuclear reactions are so closely connected with Photo-nuclear reactions that sometimes they are often 
called "Photo-nuclear". The one-photon exchange mechanism dominates in Electro-nuclear reactions, and the 
electron can be substituted by a flux of photons. Folding this flux with the gamma-nuclear cross-section, we 
arrive at an acceptable description of the electro-nuclear physics. 

Final states in gamma and electro nuclear physics are described using chiral invariant phase-space decay at low 

gamma or equivalent photon energies, and quark gluon string model at high energies. 

We will present the modeling of this physics in Geant4, and show results from practical applications. 



1. Gamma nuclear reactions 
1.1. Cross-section calculation 

The Photo-nuclear cross sections parameterized in 
the Geant4 covers all incident photon energies from 
hadron production threshold up-wards. The parame- 
terization is subdivided into five energy regions, each 
corresponding to the physical process that dominates. 
These are the Giant Dipole Resonance (GDR) region, 
the "quasi-deuteron" region, the A isobar region char- 
acterized by the dominant peak in the cross section 
which extends from the pion threshold to 450 MeV, 
the Roper resonance region that extends from roughly 
450 MeV to 1.2 GeV, and the Reggeon-Pomeron re- 
gion above 1.2 GeV. 

From the Geant4|2^| Photo- nuclear data base 
currently 14 nuclei are used in defining the parame- 
terization: 1 H, 2 H, 4 He, 6 Li, 7 Li, 9 Be, 12 C, 16 0, 27 Al, 
40 Ca, Cu, Sn, Pb, and U. The result is a function 
of A and e = log{E~), where E 7 is the incident 
photon energy. The cross-section is the sum of the 
components which parameterize each energy region. 

The cross section in the GDR region can be de- 
scribed as the sum of two peaks, 

GDR(e) = th(e, &i, si) • exp{c\ — p\ ■ e) + 
th(e, b 2 , s 2 ) ■ exp(c 2 — p% ■ e). 

The exponential describes the falling edge of the res- 
onance which has power law behavior. The function 

th(e, bi, Si) = — 1 (1) 

1 + exp{^f) 



describes the rising edge of the resonance. It is the 
nuclear-barrier-reflection function and behaves like a 
threshold, cutting off the exponential. The exponen- 
tial powers pi and p 2 are 

pi = l,p 2 = 2 for A < 4 
pi = 2,p 2 = 4 for 4 < A < 8 
Pl = 3,p 2 = 6 for 8 < A < 12 
Pi = 4,J3 2 = 8 for A > 12. 

The ^4-dependent parameters bi, Ci and Si were found 
for each of the 14 nuclei and are interpolated for 
other nuclei. 

The A isobar region was parameterized as 

* / j f s d-th(e,f,g) 

A(e,d,f,g,r,q) = — ^, (2) 

1 + r ■ (e — qY 

where d is an overall normalization factor, q can be 
interpreted as the energy of the A isobar and r as the 
inverse A width, th is the threshold function. The 
^4-dependence of these parameters is as follows: 

• d = 0.41 • A (for !H it is 0.55, for 2 H it is 0.88), 
i.e. the A yield is proportional to A; 

• / = 5.13 — .00075 • A. exp(f) shows how the 
pion threshold depends on A. 

• g = 0.09 for A > 7 and 0.04 for A < 7; 

• q = 5.84 — 1 QQ 3 . A i , which means that the 
"mass" of the A isobar moves to lower energies; 

• r = 11.9 - 1.24 • log(A). r is 18.0 for 1 H. 
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The A-dependence of the /, q and r parameters is 
due to the A + N — > N + N reaction, which can take 
place in the nuclear medium below the pion threshold. 

The quasi-deuteron contribution was parameterized 
with the same form as the A contribution but without 
the threshold function: 



QD{e, v, w, u) 



1 + w • (e — u) 2 



(3) 



For 1 H and 2 H the quasi-deuteron contribution is al- 
most zero. For these nuclei the third baryonic reso- 
nance was used instead, so the parameters for these 
two nuclei are quite different, but trivial. The param- 
eter values are given below. 



• v 



ezp(-1.7+n-0.84) 
l+eap(7-(2.38-o)) ' 



where a — log (A). The 
A-dependence in the quasi-deuteron region is 
stronger than A°- s4 It contributes only little for 
light nuclei. For 1 H it is 0.078 and for 2 H it is 
0.08. 

• u = 3.7 and w — 0.4. Experimental informa- 
tion is insufficient to fix their ^-dependence. For 
both 1 H and 2 H we have u = 6.93 and w = 90. 

The Roper contribution was parameterized using the 
same form as the quasi-deuteron contribution: 



Tr(e, v, to, u) 



(4) 



1 + w ■ (e — u) 2 
Using a = log(A), the values of the parameters are 

• v = exp(-2. + a-0.84). For *H it is 0.22 and for 
2 H it is 0.34. 

• u = 6.46 + a- 0.061 (for X H and for 2 H it is 6.57), 
so the Roper mass increases with A. 

• w = 0.1 + a ■ 1.65. For X H it is 20.0 and for 2 H 
it is 15.0. 

The Regge-Pomeron contribution was parametrized in 
terms of two exponentials describing the Pomeron and 
higher Reggeon contributions respectively: 

RP(e,h) = /i-i/i(7.,0.2)-(0.0116-e( e016 )+0.4- e (- e0 - 2 )), 

(5) 

with h = A- exp(-a ■ (0.885 + 0.0048 • a)). 



is applied to quarks, and simple kinematic mecha- 
nisms are used to model the hadronization of quarks 
into hadrons. Along with relativistic kinematics and 
the conservation of quantum numbers, the following 
concepts are introduced: 

The Quasmon is any excited hadronic system, and 
can be viewed as a generalized hadron. At the con- 
stituent level, a Quasmon may be thought of as a 
bubble of quark-parton plasma in which the quarks 
are massless. The quark-partons in the Quasmon are 
massless and homogeneously distributed over the in- 
variant phase space. It may also be considered as a 
bubble of the three-dimensional Feynman- Wilson 0] 
parton gas. The traditional hadron is a particle de- 
fined by quantum numbers and a fixed mass or width. 
The quark content of the hadron is a secondary con- 
cept constrained by the quantum numbers. The Quas- 
mon, however, is defined by its quark content and 
mass, and the concept of a well defined particle with 
quantum numbers is of secondary importance. 

The quark fusion hypothesis determines the rules 
of final state hadron production, with energy spectra 
reflecting the momentum distribution of the quarks in 
the system. Fusion occurs when a quark-partons in a 
Quasmon join to form a hadron. In cases of multiple 
Quasmon, quark-partons may be exchanged between 
the two Quasmons. Resulting hadrons are constrained 
to be produced on mass shell. The type of the out- 
going hadron is selected using combinatoric and kine- 
matic factors consistent with conservation laws. 

The only non-kinematic concept is the hypothesis of 
a critical temperature of the Quasmon. This has a 35- 
year history, starting with Ref. 0] and is based on the 
experimental observation of regularities in the inclu- 
sive spectra of hadrons produced in different reactions 
at high energies. Qualitatively, the critical tempera- 
ture hypothesis assumes that the Quasmon cannot be 
heated above a certain temperature. Adding more en- 
ergy to the system increases only the number of con- 
stituent quark-partons while the temperature remains 
constant. The critical temperature T c = 180 — 200 
MeV is the principal parameter of the model. The 
choice of this parameter is motivated from the results 
shown in Fig [31 

For the sake of briefness, we will only include the 
solution of the vacuum problem in this paper, and 
refer for the solution of the in-medium equations to 
the CHIPS publications [13, [13 • 



1.2. Final state generation, low energies 

For simulating final states for gamma-nuclear reac- 
tions, we are using the chiral invariant phase-space 
(CHIPS) approach]g[l(j[ll|. 

The CHIPS model uses a set of simple rules which 
govern microscopic quark-level behavior to model 
macroscopic hadronic systems. The invariant phase 
space distribution as a paradigm of thcrmalizcd chaos 



1 .2.1 . Solution of the vacuum equations. 

To generate hadron spectra from free Quasmons, 
the number of partons in the system must be found. 
For a finite system of N partons with a total invari- 
ant mass M, the invariant phase space integral, <I>Ar, is 
proportional to M 2N ~ 4 . At a temperature T the sta- 

Af 

tistical density of states is proportional to e~~ and 
the probability to find a system of N quark-partons 
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in a state with mass M is dW oc M 2N 4 e t dM. For 
this kind of probability distribution the mean value of 
M 2 is 

< M 2 >= 4TV(TV — 1) • T 2 . (6) 

For large TV we obtain for massless particles the well- 
known < M >= V< M 2 > = 2NT result. 

After a nucleon absorbs an incident quark-parton, 
such as a real or virtual photon, for example, the 
newly formed Quasmon has a total of TV quark- 
partons, where TV is determined by Eq. Choosing 
one of these quark-partons with energy k in the center 
of mass system (CMS) of TV partons, the spectrum of 
the remaining TV — 1 quark-partons is given by 



kdk " (M ^' 



2JV— 6 



(7) 



where Mn_i is the effective mass of the N — 1 quark- 
partons. The effective mass is a function of the total 
mass M. 



M 2 - 2kM, (8) 

so that the resulting equation for the quark-parton 
spectrum is: 



dW (x{1 _2_k )N _ 3 



(9) 



kdk v M' 
In order to decompose a Quasmon into a hadron and 
a residual Quasmon, one needs the probability of two 
quark-partons to produce the effective mass of the 
hadron. We calculate the spectrum of the second 
quark-parton by following the same argument used 
to determine Eq. |5J One quark-parton is chosen from 
the residual TV — 1. It has an energy q in the CMS of 
the TV — 1 quark-partons. The spectrum is obtained 
by substituting TV — 1 for TV and M N -\ for M in Eq.0 
and then using Eq. |5| to get 

N-4 



dW 
qdq 



1 - 



2<? 



Mjl 



■2k_ 

M 



(10) 



To ensure that the fusion will result in a hadron of 
mass fj,, we apply the mass shell constraint for the 
outgoing hadron, 

/i 2 = 2— F Ji= -q- (l-cos0). (11) 

Here 9 is the angle between the momenta k and q of 
the two quark-partons in the CMS of TV — 1 quarks. 
The kinematic quark fusion probability for any pri- 
mary quark-parton with energy k is then: 



N-4 



2q 



x 8 



2kq(l 



Using the (5-function to perform the integration over 
q one gets: 



P{k,M,tx) 



Mk{l-cos6) 

2 



N-4, 



,,2 /i 2k 



2k(l — cos 



M -2k 
x d I 1- 



P(k,M,n) 



After the substitution z 
this becomes 



1 



Mk(l - cost 

2 



Mk{l-cos9) 
= 1-^2_ = 1 



N-4 



(14) 



Mjv-i 



Mk(l-cos6) ' 



P(fc,M, M ) 



M -2k 
4fc 



(15) 



where the limits of integration are when cos 9 
1 -JTk^ and 



= 1 - 



(16) 



when cos 6 = — 1. The resulting range of 9 is therefore 

2 

— 1 < cos 9 < 1 — . Integrating from to z yields 

M - 2k 



,N-3 



4k- (N — 3) 



(17) 



and integrating from to z max yields the total kine- 
matic probability for hadronization of a quark-parton 
with energy k into a hadron with mass [i: 



M — 2k 
4k ■ (TV - 3) 



z, 



N-3 



(18) 



The ratio of expressions El an d El can be treated as 
a random number, R, uniformly distributed on the 
interval [0,1]. Solving for z then gives 



z = v R ■ z n 



(19) 



In addition to the kinematic selection of the two 
quark-partons in the fusion process, the quark con- 
tent of the Quasmon and the spin of the candidate fi- 
nal hadron are used to determine the probability that 
a given type of hadron is produced. Because only the 
relative hadron formation probabilities are necessary, 
overall normalization factors can be dropped. Hence 
the relative probability can be written as 



P h (k,M,v) = (2s h + l)-z^- x 3 -C%. 



(20) 



1 _ 2fe 
1 At 



qdqd cos 0(12) Here, only the factor z^ ax 3 is used since the other fac- 
tors in eauationll8larc constant for all candidates for 
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the outgoing hadron. The factor 2sh + 1 counts the 
spin states of a candidate hadron of spin s/ l; and Cq 
is the number of ways the candidate hadron can be 
formed from combinations of the quarks within the 
Quasmon. In making these combinations, the stan- 
dard quark wave functions for pions and kaons were 
used. For n and n' mesons the quark wave functions 
T] = - ^ and i = 5H±^ + J| were used. No 

mixing was assumed for the u> and <j> meson states, 
hence u = uu +^ d and ip = ss. 



1.3. Final state generation, high energies 

At high energies we use quark-gluon string model 
and a diffractive Ansatz for string excitation to de- 
scribe the interactions of real and virtual photons 
with nuclei. A description of the means of doing 
this can be found in a separate paper in the present 
proceedings^. 



2. Electro-nuclear scattering 

Electro-nuclear reactions are very connected with 
Photo-nuclear reactions. They are sometimes called 
"Photo-nuclear" because the one-photon exchange 
mechanism dominates the reaction. In this sense elec- 
trons can be replaced by a flux of equivalent photons. 
This is not completely true, because at high energies 
diffractive mechanisms are possible, but these types 
of reactions are beyond the scope of this discussion. 



2.1. Electro-nuclear cross-sections 

The Equivalent Photon Approximation (EPA) was 
proposed by E. Fermi [ll and developed by C. Weiz- 
sacker and E. Williams 31 an d by L. Landau and E. 
Lifshitz The covariant form of the EPA method 
was developed in Refs. a 4i] and 5]. When using this 
method it is necessary to take into account that real 
photons are always transversely polarized while vir- 
tual photons may be longitudinally polarized. In gen- 
eral the differential cross section of the Electro-nuclear 
interaction can be written as 



d 2 a 



dydQ 2 irQ 2 
where 

1 ~y i 1 it 



{S T L-{vT + o- L )-S L -a L ), (21) 



Stl = y- 



V 2 + ^r 

L ~ 2^ ~ ~Q T '' 



The differential cross section of the Electro-nuclear 
scattering can be rewritten as 



d 2 a e 



ay 



dydQ 2 ttQ 2 



f 



y- 



91 4 

E 2 



Q 2 



o-^a, (22) 



where ct 7 »a = a 1 A(y) for small Q 2 and is written as 
a function of e, v, and Q 2 for large Q 2 . Interactions 
of longitudinal photons are normally included in the 
effective a 7 *A cross section through the e factor, 
but in the present method, the cross section of 
virtual photons is considered to be e- independent. 
The Electro-nuclear problem, with respect to the 
interaction of virtual photons with nuclei, can thus 
be split in two. At small Q 2 it is possible to use the 
<J 1 (v) cross section. In the Q 2 >> m 2 region it is 
necessary to calculate the effective ct 7 « (e, v, Q 2 ) cross 
section. 

Following the EPA notation, the differential cross 
section of Electro-nuclear scattering can be related to 
the number of equivalent photons dn = -29-, For y « 

1 and Q 2 < Am 2 the canonical method § leads to 



ydn(y) 
dy 



2a 



ln(y). 



(23) 



In 0, integration over Q 2 for v 2 » Q 2 max ~ m 2 
leads to 



ydn(y) _ (l + {l-y) 
dy 



-*„(-*_) + (!-„) 

i -y 



(24) 

In the y « 1 limit this formula converges to Ea. (|23|l . 
But the correspondence with Ea. (|23|) can be made 
more explicit if the exact integral 

ydn(y) a ( \ + {l-yf (2 - yf , 

~^y- = « { 2 h V)h ~^— h ' " 



where l\ 



In 



,„ f v 2 +QlaJE 2 \ n 2 

\V 2 +Q' 2 mzn /E 2 ) ' Wm 



(25) 



is calculated for 



Q 



4m; 



rnax(m e ) 1— y ' 

The factor (1 — y) is used arbitrarily to keep 
Qmax(m c ) > <3mm' which can be considered as a 
boundary between the low and high Q 2 regions. The 
transverse photon flux can be calculated as an integral 
of Eq. I|25|l with the maximum possible upper limit 



Q 2 max{max) = 4£ 2 (1 - y). 



(26) 



It can be approximated by 

ydnjy) = _2a ({2- y) 2 + y 2 
dy 7i 



ln( 7 ) - 1 , (27) 



where 7 = — . It must be pointed out that neither 
this approximation nor Ea. (|25|l works at y ~ 1. The 
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Figure 1: Relative contribution of small Q 2 to the 
"photon flux" for (a) 1 GeV electrons and (b) 10 GeV 
electrons. In figures (c) and (d) the photon distribution 
dn(v,Q 2 ) is multiplied by the Photo-nuclear cross section 
and integrated over Q 2 . The dashed lines are integrals 
over low Q 2 (under the dashed line in the first two 
figures), the solid lines are integrals over high-Q 2 (above 
the dashed lines in the first two figures). 



formal limit of the method is y < 1 — ■ 2 —. 

Fig. ^a,b) shows the energy distribution for the 
equivalent photons. The low-Q 2 flux is calculated us- 
ing Eq. (dashed lines) and Eq. (dotted lines) . 
The total flux is calculated using Ea. (|27Jl (the solid 
lines) and using Ea. I)25|l with the upper limit de- 
fined by Eq. (|26[1 (dash-dotted lines visible only around 
v E e ). We find that to calculate the number of low- 
Q 2 equivalent photons or the total number of equiv- 
alent photons one can use the approximations given 
by Ea. (|25|l and Eq. ljT7|) . respectively, instead of us- 
ing Ea. (|25|l . Comparing the low-Q 2 photon flux and 
the total photon flux we find that the low-Q 2 photon 
flux is about half of the the total. From the interac- 
tion point of view the decrease of ct 7 * with increasing 
Q 2 must be considered. The cross section reduction 
for the virtual photons with large Q 2 is governed by 
two factors. The cross section drops with Q 2 as the 
squared dipole nucleon form-factor 



G 2 D (Q 2 ) 



Q 2 



(843 MeVf 



(28) 



and the thresholds the reactions are shifted to 
higher v by a factor which is the difference be- 
tween the K and v values. Following the method pro- 
posed in @ , (7 7 * at large Q 2 can be approximated as 



(1 - x)a 1 {K)G 2 D {Q 2 )e b{t ' K) - r+c{t ' K) - r \ (29) 



where r 



The e-dependence of the 



i{e, K) and b(e, K) functions is weak, so for simplic- 



ity the b(K) and c(K) functions are averaged over e. 
They can be approximated as 



b(K) 



K 



185 MeV 



0.85 



,c(K) 



K 



1390 MeV 



(30) 

The integrated photon flux folded with the cross 
section approximated by Ea. l|29|) is shown in 
Fig- ^c,d). We show separately the low-Q 2 region 



{Q < 

region (Q 2 > Q 



) 2 

' max(m e 



dashed lines), and the high-Q 
solid lines). These functions 



max(m e ) 

will be integrated over ln{u), hence because of the 
Giant Dipole Resonance contribution, the low-Q 2 
part covers more than half the total eA — > hadrons 
cross section. But at v > 200 MeV , where the hadron 
multiplicity increases, the large Q 2 part dominates. 
In this sense, for a better simulation of the production 
of hadrons by electrons, it is necessary to simulate 
the high-Q 2 part as well as the low-Q 2 part. 

Taking into account the contribution of high-Q 2 
photons it is possible to use Ea. (|27|l with the over- 
estimated (T^-A = &ia{v) cross section. The slightly 
over-estimated Electro-nuclear cross section is 



'eA 



where 



(2Zn(7) - 1) • Jx ~ 



2J CM, 



Ji(E e 



J2(E e ) = - 
IT 



o- 1 A( v )dl n { v ) (32) 



J 3 (E e 



va^Aiy^dlniy), and (33) 



v 2 a 1 A{v)dln(v). (34) 



The equivalent photon energy v = yE can be obtained 
for a particular random number R from the equation 



R = 



{2ln{ 1 )-l)J 1 { V )- l ^{2J 2 { V )-^) 



(2Zn( 7 ) - l)Jx{E e ) - l -^(2J 2 (E e ) - ^) ' 



(35) 

Ea. l|25|l is too complicated for the randomization of 
Q 2 but there is an easily randomized formula which 
approximates Eq. (j25|) above the hadronic threshold 
(E > 10 MeV). It reads 



Q 



ydn(y,Q 2 ) 
dydQ 2 



dQ 2 



aD(y) 
where 

L(y,Q 2 )=ln(F(y) + (e p ^ -1 + 



L(y,Q 2 )-U(y), 



Q 2 



i 2 . ' 

rn t n 



(36) 



(37) 
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Figure 2: Integrals of Q 2 spectra of virtual photons for 
three energies 10 MeV, 100 MeV, and 1 GeV at 
y = 0.001, y = 0.5, and y = 0.95. The solid line 
corresponds to Eg. 12511 and the dashed line (which 
almost everywhere coincides with the solid line) 
corresponds to Eg. 12511 . 



D(y) = l-y + 



2 ' 



and 



U(y) = P(y) 



1 - 



Q 



2 

rain 



) 2 

* m a.r 



, with 



- max 



P(y) 



D(yY 



(38) 



(39) 



(40) 



(41) 



The Q 2 value can then be calculated as 

Q 2 



) 2 . 

■> in i n 



l-e 



P(v). 



_^ e RL(y,C 



C )-(1-R)-U(y) 



F{y) 



(42) 

where R is a random number. In Fig. [21 Eq.J^SJ (solid 
curve) is compared to Eo . 1)36(1 (dashed curve). Be- 
cause the two curves are almost indistinguishable in 
the figure, this can be used as an illustration of the 
Q 2 spectrum of virtual photons, which is the deriva- 
tive of these curves. An alternative approach is to use 
Eq. (|25|l for the randomization with a three dimen- 
sional table ^(Q 2 ,y,E e ). 

After the v and Q 2 values have been found, the value 
of <r 7 .yi(i/, Q 2 ) is calculated using Eq.J2jJ). Note that 
if R ■ a^A.{y) > o~-f*A(y, Q 2 ), no interaction occurs and 
the electron keeps going. 

2.2. Final state generation 

Final states are generated using the single photon 
exchange assumption. Sampling the equivalent pho- 



ton distribution described in the previous section al- 
lows to construct an exchange particle that in turn 
can be treated by the mechanisms used for gamma 
nuclear scattering. The question to be answered is 
that of the absorption mechanisms of these particles 
in the context of parton exchange diagrams. 

In the example of the Photo-nuclear reaction dis- 
cussed in the comparison section, namely the descrip- 
tion of 90° proton and deuteron spectra in A(j,X) 
reactions at E 1 = 59 — 65 MeV, the assumption on the 
initial Quasmon excitation mechanism was the same. 
The description of the 90° data was satisfactory, but 
the generated data showed very little angular depen- 
dence, as the velocity of Quasmons produced in the 
initial state was small, and the fragmentation process 
was almost isotropic. Experimentally, the angular de- 
pendence of secondary protons in photo-nuclear re- 
actions is quite strong even at low energies (see, for 
example, Ref. [25 )■ This is a challenging experimen- 
tal fact which is difficult to explain in any model. It's 
enough to say that if the angular dependence of sec- 
ondary protons in the 7 40 Ca interaction at 60 MeV 
is analyzed in terms of relativistic boost, then the 
velocity of the source should reach 0.33c; hence the 
mass of the source should be less than pion mass. The 
main subject of the present publication is to show that 
the quark-exchange mechanism used in the CHIPS 
model can not only model the clusterization of nu- 
cleons in nuclei and hadronization of intra-nuclear ex- 
citations into nuclear fragments, but can also model 
complicated mechanisms of interaction of photons and 
hadrons in nuclear matter. 

Quark-exchange diagrams help to keep track of the 
kinematics of the quark-exchange process To apply 
the mechanism to the first interaction of a photon 
with a nucleus, it is necessary to assume that the 
quark-exchange process takes place in nuclei contin- 
uously, even without any external interaction. Nu- 
cleons with high momenta do not leave the nucleus 
because of the lack of excess energy. The hypothesis 
of the CHIPS model is that the quark-exchange forces 
between nucleons 15] continuously create clusters in 
normal nuclei. Since a low-energy photon (below the 
pion production threshold) cannot be absorbed by a 
free nucleon, other absorption mechanisms involving 
more than one nucleon have to be used. 

The simplest scenario is photon absorption by a 
quark-parton in the nucleon. At low energies and in 
vacuum this does not work because there is no corre- 
sponding excited baryonic state. But in nuclear mat- 
ter there is a possibility to exchange this quark with a 
neighboring nucleon or a nuclear cluster. The diagram 
for the process is shown in the upper part of Fig. 0] 
In this case the photon is absorbed by a quark-parton 
from the parent cluster PCi, and then the secondary 
nucleon or cluster PC2 absorbs the entire momentum 
of the quark and photon. The exchange quark-parton 
q restores the balance of color, producing the final- 
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Figure 3: Average meson multiplicities in 
proton-antiproton and in electron-positron annihilation, 
as a function of the CMS energy of the interacting 
hadronic system. The points are experimental data [73| 
and the lines are CHIPS Monte Carlo calculations at 
several values of the critical temperature parameter T. 



state hadron F and the residual Quasmon RQ. The 
process looks like a knockout of a quasi-free nucleon 
or cluster out of the nucleus. It should be emphasized 
that in this scenario the CHIPS event generator pro- 
duces not only "quasi-free" nuclcons but "quasi-free" 
fragments too. The yield of these quasi-free nucleons 
or fragments is concentrated in the forward direction. 

The second scenario, shown in the lower part of 
Fig0]which provides for an angular dependence is the 
absorption of the photon by a colored fragment (CF2 
in Fig. |3J) ■ In this scenario, both the primary quark- 
parton with momentum k and the photon with mo- 
mentum q 7 are absorbed by a parent cluster (PC2 
in the lower part of Fig. 0J, and the recoil quark- 
parton with momentum q cannot fully compensate the 
momentum k + q~. As a result the radiation of the 
secondary fragment in the forward direction becomes 
more probable. 

In both cases the angular dependence is defined by 
the first act of hadronization. The further fragmenta- 
tion of the residual Quasmon is almost isotropic. 

It was shown in above that the energy spectrum of 
quark partons in a Quasmon can be calculated as 



dW 

k*dk* 



oc 



f 



2k* 
~M 



N-3 



(43) 



where k* is the energy of the primary quark-parton in 
the Center of Mass System (CMS) of the Quasmon, 
M is the mass of the Quasmon, and N, the number 
of quark-partons in the Quasmon, can be calculated 
from the equation 

< M 2 >= 4 • N ■ (N - 1) ■ T 2 . (44) 

Here T is the temperature of the system. 




1A 
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PC 2 (A 2 ) CP 



pc #iL 
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CF 



>- Q. ^ > O F > 




Figure 4: Diagrams of photon absorption in the quark 
exchange mechanism. PCi,a are the parent clusters with 
bound masses involved in the quark-exchange. 
CFi,2 are the colored nuclear fragments. F(/i) is the 
outgoing hadron with mass fi. RQ is the residual 
Quasmon which carries the rest of the excitation energy 
and momentum. M m i n characterizes its minimum mass 
defined by its quark content. Dashed lines indicate 
colored objects. The photon is absorbed by a 
quark-parton k from the parent cluster PCi. 



In the first scenario of the jA interaction (Fig. 0J, 
as both interacting particles are massless, we assumed 
that the cross section for the interaction of the pho- 
ton with a particular quark-parton is proportional to 
the charge of the quark-parton squared, and inversely 
proportional to the mass of the photon-parton system 
s, which can be calculated as 



2wfc(l - cos(0 fc )). 



(45) 



Here u) is the energy of the photon, and k is the energy 
of the quark-parton in the Laboratory System (LS): 



En +Pn ■ cos( 
Mn 



(46) 



In the case of a virtual photon, equation l|45|) can be 
written as 



2k(u> - q 1 ■ cos(0 k )), 



(47) 



where g 7 is the momentum of the virtual photon. In 
both cases equation (|43|l transforms into 



dW 
dk 1 



oc 



1 



2k* 



N-3 



(48) 



and the angular distribution in cos(0fc) converges to a 
J-function: in the case of a real photon cos^) = 1, 
and in the case of a virtual photon cos(0k) = 
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Figure 5: Comparison of the CHIPS model results (lines 
in the figure) with the experimental data [rs| on proton 
spectra at 90° in the Photo- nuclear reactions on 40 Ca at 
59-65 MeV (open circles), and proton spectra at 60° 
(triangles) and 150° (diamonds). Statistical errors in the 
CHIPS results are not shown but can be judged by the 
point-to-point variations in the lines. The comparison is 
absolute, using the value of total Photo- nuclear cross 
section of 5.4 mb for Ca, as given in Ref. |Tgfl . 



In the second scenario for the photon interaction 
(lower part of Fig.^J we assumed that both the photon 
and the primary quark-parton, randomized according 
to equation l|43[) . enter the parent cluster PC2, and af- 
ter that the normal procedure of quark exchange con- 
tinues, in which the recoiling quark-parton q returns 
to the first cluster. 

An additional parameter in the model is the rela- 
tive contribution of both mechanisms. As a first ap- 
proximation we assumed equal probability, but in the 
future, when more detailed data are obtained, this pa- 
rameter can be adjusted. 



3. Comparison with experiment 

We begin the comparison with the data on proton 
production in the 40 Ca(7, X) reaction at 90° at 59-65 
MeV [H, and at 60° and 150° at 60 MeV 0]. We 
analyzed these data together to compare the angular 
dependence generated by CHIPS with experimental 
data. The data are presented as a function of the in- 



Figure 6: Comparison of the CHIPS model results (lines 
in the figure) with the experimental data on proton 
spectra at 57°, 77°, 97°, 117°, and 127° in the 
Photo-nuclear reactions on 12 C at 123 MeV (open 
circles). The value of the total Photo- nuclear cross 
section was set at 1.8 mb. 
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C(y,p) reaction at E = 1 51 MeV 




variant inclusive cross section / = 



da 
p p dE p 
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Figure 7: Same as in Fig. |H| for the photon energy 151 
MeV. 
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Figure 8: Comparison of the CHIPS model results (line 
in the figure) with the experimental data [23| (open 
circles) on the proton spectrum measured in parallel 
kinematics in the 12 C(e,e'p) reaction at an energy 
transfer equal to 210 MeV and momentum transfer equal 
to 585 MeV/c. Statistical errors in the CHIPS result are 
not shown but can be judged by the point-to-point 
variations in the line. The relative normalization is 
arbitrary. 



on the variable k — Tp + Pp , 



where T p and p p are the 



kinetic energy and the momentum of the secondary 
proton. As one can see from Fig. [5] the angular de- 
pendence of the proton yield in photo-production on 
40 Ca at 60 MeV is reproduced quite well by the CHIPS 
event generator. 

The second set of measurements that we use for the 
benchmark comparison deals with the secondary pro- 
ton yields in 12 C(7, X) reactions at 123 and 151 MeV 
|22|. which is still below the pion production thresh- 
old on a free nucleon. Inclusive spectra of protons 
have been measured in 7 12 C reactions at 57°, 77°, 
97°, 117°, and 127°. Originally, these data were pre- 
sented as a function of the missing energy. We present 
the data in Figs. El and [7| together with CHIPS cal- 
culations in the form of the invariant inclusive cross 
section dependent on k. 

The agreement between the experimental data and 
the CHIPS model results is quite remarkable. Both 
data and calculations show significant strength in the 
proton yield cross section up to the kinematic limits 
of the reaction. The angular distribution in the model 
is not as prominent as in the experimental data, but 



agrees well qualitatively. 

Using the same parameters, we applied the CHIPS 
event gen erator to the 12 C(e,e'p) reaction measured 
in Ref . [23 ■ The proton spectra were measured in par- 
allel kinematics in the interaction of virtual photons 
with energy w = 210 MeV and momentum g 7 = 585 
MeV/c. To account for the experimental conditions in 
the CHIPS event generator, we have selected protons 
generated in the forward direction with respect to the 
direction of the virtual photon, with the relative angle 
<d qp < 6°. The CHIPS generated distribution and the 
experimental data are shown in Fig. [H] in the form of 
the invariant inclusive cross section as a function of k. 
The CHIPS event generator works only with ground 
states of nuclei so we did not expect any narrow peaks 
for 1 p 3 / 2 -shell knockout or for other shells. Neverthe- 
less we found that the CHIPS event generator fills 
in the so-called " 1 Si/2-shell knockout" region, which 
is usually artificially smeared by a Lorentzian |24|. 
In the regular fragmentation scenario the spectrum 
of protons below k = 300 MeV is normal; it falls 
down to the kinematic limit. The additional yield at 
k > 300 MeV is a reflection of the specific first act 
of hadronization with the quark exchange kinematics. 
The slope increase with momentum is approximated 
well by the model, but it is obvious that the yield 
close to the kinematic limit of the 2^2 reaction can 
only be described in detail if the excited states of the 
residual nucleus are taken into account. 

The angular dependence of the proton yield in low- 
energy photo-nuclear reactions is described in the 
CHIPS model and event generator. The most impor- 
tant assumption in the description is the hypothesis of 
a direct interaction of the photon with an asymptot- 
ically free quark in the nucleus, even at low energies. 
This means that asymptotic freedom of QCD and dis- 
persion sum rules | l7| can in some way be generalized 
for low energies. The knockout of a proton from a 
nuclear shell or the homogeneous distributions of nu- 
clear evaporation cannot explain significant angular 
dependences at low energies. 

The same mechanism appears to be capable of mod- 
eling proton yields in such reactions as the 16 C(e,e'p) 
reaction measured at MIT Bates [2^, where it was 
shown that the region of missing energy above 50 MeV 
reflects "two-or-more-particle knockout" (or the "con- 
tinuum" in terms of the shell model). The CHIPS 
model may help to understand and model such phe- 
nomena. 
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